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Quantum Metamagnetic Transitions Induced by Changes in Fermi-Surface Topology 
—Applications to a Weak Itinerant-Electron Ferromagnet; ZrZn 2 
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We clarify that metamagnetic transitions in three dimensions show unusual properties as 
quantum phase transitions if they are accompanied by changes in Fermi-surface topology. An 
unconventional universality deeply affected by the topological nature of Lifshitz-type transi- 
tions emerges around the marginal quantum critical point (MQCP). Here the MQCP is defined 
by the meeting point of the finite temperature critical line and a quantum critical line running 
on the zero temperature plane. The MQCP offers a marked contrast with the Ising universal- 
ity and the gas-liquid-type criticality satisfied for conventional metamagnetic transitions. At 
the MQCP, the inverse magnetic susceptibility \ X nas diverging slope as a function of the 
magnetization m (namely, \alx~ 1 /dm\ — > oo) in one side of the transition, which should not 
occur in any conventional quantum critical phenomena. The exponent of the divergence can 
be estimated even at finite temperatures. We propose that such an unconventional universality 
indeed accounts for the metamagnetic transition in ZrZn2. 
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Itinerant ferromagnets such as ZrZn.2, 1 ' 2 UGe2 3 and 
nearly ferromagnetic metals such as Sr3Ru207 show 
metamagnetic transitions. The magnetizations show 
jumps at magnetic fields separating the low-field phase 
with a smaller magnetic moment from the high-field 
phase with a higher moment. The first-order transition 
characterized by the magnetization jump terminates at a 
finite-temperature critical point. The universality around 
the critical point is regarded as the Ising type, which 
is equivalent to the gas-liquid critical points. The criti- 
cal temperature can, however, be controlled to zero, for 
example, by pressure, which offers a field of quantum 
critical phenomena. The metamagnetic transition and its 
quantum critical point (QCP) in itinerant electron sys- 
tems have attracted interest because of its intriguing na- 
ture of fluctuations leading to possible non- Fermi-liquid 
behavior as well as superconductivity found in UGe2 near 
the metamagnetic transition. 5 

A simple scope of the QCP is offered by suppres- 
sions of critical temperatures of gas-liquid (or equiva- 
lently Ising) transitions, T c , by quantum fluctuations. 6 
The QCP of magnetic transitions in itinerant electron 
systems has been analyzed by the same framework of 
the suppressed symmetry breaking by the quantum fluc- 
tuations, 7-9 which is expressed essentially by the d + z 
dimensional Ising criticality with the spatial dimension- 
ality d being added by the dynamical exponent z repre- 
senting the quantum dynamics. The quantum criticality 
of the metamagnetic transition has been interpreted so 
far by the same conventional QCP of the Ising type. 10 

Recently, however, for two-dimensional systems, a 
completely different type of quantum critical phenom- 
ena has been proposed for the metal-insulator transi- 
tions 11 ' 12 as well as for the Lifshitz transition, 13 where 
the topological nature of the transition deeply modifies 
the above conventional picture. In this case, the phase 
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Fig. 1. (a) Schematic phase diagram around the marginal quan- 
tum critical point (MQCP), where a quantum critical line (QCL) 
and a finite-temperature critical line at T c meet. Here, B controls 
quantum fluctuations and corresponds to, for instance, magnetic 
field. Experimentally, pressure P may control the ratio of itiner- 
ancy (bandwidth) W to the interaction U. (b) Schematic phase 
diagram proposed for ZrZn2. x An extended range of applied pres- 
sure is illustrated including the negative side of P, where an 
additional metamagnetic transition has been speculated. TCP 
represents the tricritical point, where the critical line at B = 
terminates and the transition becomes first order at B = below 
the tricritical temperature. This first-order line below TCP is the 
edge of the metamagnetic "wing" , which is the main subject of 
this letter. The dashed lines were regarded just as a crossover in 
Ref. 1 based on the picture of the conventional QCP. 



diagram is qualitatively different from that of the usual 
QCP, as is shown in Fig. 1(a). Although the critical line 
at nonzero critical temperature T c is characterized by 
the conventional Ising universality, there emerges a quan- 
tum critical line (QCL) of topological transitions starting 
from the QCP on the zero-temperature plane. This QCP 
is called marginal quantum critical point (MQCP) as the 
zero-temperature critical line terminates at this point. 

In this letter, we propose that the metamagnetic tran- 
sition in three dimensional systems also bears an un- 
conventional feature characterized neither by the Ising 
universality nor by any other universality of symmetry- 
breaking transitions, when the metamagnetic transition 
is involved and triggered by the topology change of the 



*E-mail:yamaji@solis.t.u-tokyo.ac.ip 



2 J. Phys. Soc. Jpn. 

Fermi surface (ex. collapse of a neck of the Fermi surface 
or emergence/disappearance of a Fermi pocket, namely 
the Lifshitz transitions 14 ). This class of metamagnetic 
critical point is a manifestation of the MQCP. We pro- 
pose that the metamagnetic transition of ZrZn2 indeed 
belongs to this unconventional class. The structure of the 
phase diagram proposed for ZrZn2 1,2 shown in Fig. 1(b) 
is consistent with our interpretation shown in Fig. 1(a) 
for the part of the metamagnetic transition beyond the 
tricritical point (TCP). Furthermore, we propose that 
the dashed lines in Fig. 1(b) are not crossover lines as 
in the interpretation in Refs.l, and 2 but represent real 
phase transition boundaries (QCL) of Lifshitz- type topo- 
logical transitions. We show that the magnetic field de- 
pendence of the inverse magnetic susceptibility x _1 is a 
key quantity for the identification of the universality of 
the MQCP, where the slope diverges toward the MQCP. 
X^ 1 observed for ZrZn2 by Uhlarz et al. 2 indeed shows 
a remarkable enhancement. Even at finite temperatures, 
we show that the MQCP universality is easily observable 
beyond a crossover from the Ising one. 

We first discuss the itinerant electron metamagnetism 
by using a mean-field theory after Wohlfarth, Rhodes, 15 
and Shimizu 16 (WRS). Through a mean-field decoupling 
of the Hubbard hamiltonian in the standard notation 

« = E ( £ ( fc ) - ^ 9J y b ) 4.4. + E m ^ « 

her i 

where n is +1 (-1) for the spin a =\ (J.), mean-field band 
dispersion for partially polarized quasiparticles, 

E k * = E(k)-r)(Um + gn B B)/2, (2) 

are obtained. The magnetization, m = nj — n^ is self- 
consistently determined afterwards. Here, the electron 
densities of each spin sub-bands, n^, n^, are given as 

/+oo 
p(e){l + exp((e-/i CT )/r)}- 1 de, (3) 
-oo 

where p(e) is the density of states (DOS) for the bare 
dispersion, e(k) and p a = p+n(U m+gpsB) /2. Then the 
magnetic-field dependence of the magnetization m is self- 
consistently calculated from Eqs. (2) and (3) with m = 
n-f — ni . For the canonical ensemble, which is justifiable 
under the charge neutrality, the total electron density, 
n = rif + n i is fixed. In the following part of the present 
letter, the chemical potential p is chosen to keep the 
electron density fixed. 

At zero temperature, the energy density 

f° = E - AV)P(£R + j(m 2 + n 2 ) (4) 

describes the nature of phase transitions. In the present 
theory, the metamagnetic instability is signaled by the 
divergence of the uniform magnetic susceptibility, \ — 
dm/dB = 2p e g(p) (1 — Up e g(p))~ 1 where p c g is the ef- 
fective DOS in the partially polarized state given as 

p cS (p) = Mm)p(h)/(p(m) + P(H))- The WRS theory 
is essentially based on the Ginzburg-Landau expansion of 
/o with respect to m, so-called the t^ 6 -theory. However, 
we note that the Ginzburg-Landau expansion is not jus- 
tified around the Lifshitz transition as we show below. 
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Fig. 2. (a) Fermi surface around a saddle point in the discon- 
nected phase (see text), (b) Fermi surface in the connected phase. 



Now we focus on the situation that necks of the Fermi 
surface can be disconnected or connected in 3D Brillouin 
zone by changing magnetic fields or applied pressure, as 
is expected in ZrZn2. 19 In this case, £l saddle point of the 
band dispersion, fcz, is located near the Fermi surface. 
Around the saddle point in three dimensions, we employ 
the axisymmetric expansion of the band dispersion as 

e(k) = e{k L ) + ^q\ + -p + 0(q% (5) 

where q = k — Ul and m\ • m* z < 0. The Fermi level 
at the Lifshitz transition is denoted as ul = s(kL) and 
the corresponding magnetic field as Bl- When the ef- 
fective mass for the xy-plane, m\ is positive (negative), 
we call the neck of the Fermi surface the electron (hole) 
neck. We also call the paramagnetic metallic phase with 
(p — Pl)tti* > (< 0), the disconnected (connected) 
phase (see Fig. 2). At p = pl, the density of states 
p(e) has a singularity originating from the changes in 
the Fermi-surface topology. For the connected-phase side 
(p — pi)m* < 0, the DOS p{e) may be expanded as 

p(e)=p(p L )+p^C + 0(C 2 ), (6) 

where £ = e — u^. On the other hand, for the 
disconnected-phase side, p(e) may be expanded as 

p(e) = P {p L ) + P { 1 /2) \C\ 1/2 + P W C + 0(C 3/2 ), (7) 

which is analytic with respect to ICI 1 ^ 2 except for the 
point C = 0. The coefficient p^ 2 ^ has a negative value. 

Substituting Eqs. (6) and (7) into Eq. (4), we expand 
fo around the Lifshitz-transition point with respect to 
5m = — sign(m*)(m — m^) and SB = B — Bl as, 

, , , (1 - Up eS (fj, L )) 2 9VB XU , v 

jo ~ const H — dm oB(m — mL) 

4:PcS(Pl) 2 

f c b/2 \5m\V 2 + 0(5m 3 ) (Sm < 0) , > 

+ \ c 3 Sm 3 + 0(Sm 4 ) (Sm > 0) ' 1 ' 

where c 5 / 2 and C3 are constants given by p(pi), p(pi), p± 

and p^+ 2 \ The region Sm < (Sm > 0) corresponds to 
the disconnected (connected) phase. Here, the Lifshitz 
transition occurs at Sm = 0. The nonanalytic term in Eq. 
(8) for the disconnected side, c 5 / 2 \Sm\ 5 ^ 2 , originates from 

the integral of the form J^ 6 " 1 cZee |e| 1/ ' 2 , which inevitably 
appears in f substituted with Eq. (7). 

Generally, the Lifshitz transitions may occur inde- 
pendently of possible metamagnetic transitions at zero 
temperature. However, if p e s has a maximum at the 
Lifshitz-transition point p = p^ as a function of p, 
X oc (1 — Up c e(p))~ 1 given below Eq.(4) diverges first 
at p = pl with the increase of U. Then a metamagnetic 
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transition may occur simultaneously with the Lifshitz 
transition at p — pl- If 1/ p e s{ph) < U and B ^ are 
satisfied, this simultaneous transition occurs as a first- 
order transition. The first-order transition terminates at 
a nonzero critical temperature T c , which has the conven- 
tional Ising universality. The first-order transition ter- 
minates at T = for U p e g{pL) - * L which is nothing 
but the MQCP. The Lifshitz-type QCL starts from the 
MQCP in the region Up c e(pL) < 1- The emergence of 
MQCP and first-order transitions is a consequence of the 
electron correlation effect, while the Lifshitz transitions 
on the QCL remain "transitions of the 2| order" pro- 
posed for non- interacting systems by Lifshitz. 14 

From the minimization of /q in Eq. (8), the singularity 
of the magnetization around the MQCP is obtained as 

-\B L -B\V S < (B<B L ) 
\B-B L \V*> (B>B L ) ' 



Sm oc 



(9) 



where, for the electron-necks case to* < (hole-necks 
case to* > 0), the critical exponent <5 is <5< — 3/2 (<5< — 
2) for B < B L and <5> = 2 (<5> = 3/2) for B > B L , which 
are clearly different from the mean-field Ising value, S = 
3. In the case of the electron necks, the singularity of the 
inverse susceptibility x _1 = dB /dm is given as 



X 1 <x 



\B L — i?) 1 / 3 oc IJml 1 / 2 (B<B L ) 
\B - B L \ l l 2 oc | <5m | (B>B L ) 



(10) 



whereas, in the case of the hole necks, x" 1 01 \Sm\ for 
B < B L and x -1 oc I^toI 1 / 2 for B > B L . At a con- 
ventional symmetry-breaking transition, x _1 oc l^mj 15-1 
vanishes with an exponent 8 > 3, resulting in a concave 
curve for the plot of x" 1 as a function of to. This applies 
to the present case with a nonzero T c as well. In marked 
contrast, around the MQCP, x _1 vanishes with a convex 
dependence with (5 = 3/2 or at most linearly with (5 = 2, 
which reflect the topological character of the transition. 
Along the QCL, on the other hand, x _1 does not vanish. 

We now examine whether the MQCP is relevant to 
understand the experimental results for ZrZii2. ZrZn2 is 
one of the canonical weak itinerant-electron ferromagnets 
with a small ordered moment 0.17 ps per Zr atom and 
low Curie temperature T c — 28. 5K, at ambient pressure. 
Recent experimental results show that the magnetization 
suddenly vanishes at P — 1.65GPa and exhibits a first- 
order transition at B — below the temperature of a 
TCP. 2 In addition, two successive metamagnetic behav- 
iors are detected for non-zero magnetic fields, 1 ' 2 as shown 
in Fig. 1(b), where one of the transitions is connected to 
the TCP at B = and P > 0. 

Several studies have also been performed to eluci- 
date the Fermi-surface topology of ZrZii2. Singh and 
Mazin 17 proposed the neck-free paramagnetic state, 
whereas electron-positron annihilation radiation mea- 
surements reveal the presence of the necks at ambient 
pressure. 18 A recent accurate ab initio electronic struc- 
ture calculation by Harima shows that the necks around 
the L point in the band closest to the Fermi surface 
change their topology depending on the location of the 
Fermi level. 19 These results together suggest that ZrZii2 
is in a subtle proximity of the Lifshitz transition. 
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Fig. 3. (a) DOS p(e) of ZrZn2 calculated by using a full potential 
LAPW method based on LDA with lattice parameters in a para- 
magnetic phase, at ambient pressure. 19 (b) An enlargement of 
(a) around e = fi^. The solid curve represents DOS fitted to the 
calculated result represented by the (red) cross, and is used in 
the present letter (see text) . The Fermi level in the paramagnetic 
phase, fio, is located below fiij. 



Here we assume that an applied pressure just renor- 
malizes the overall energy scale such as U/W and leave 
the band structure calculated by LDA 19 unchanged. 
Then saddle points associated with electron necks sur- 
rounding the L point of the cubic Laves phase, around 
which quadratic band dispersions such as Eq. (5) are ex- 
pected, are located above the paramagnetic Fermi level 
as is shown in Fig. 3. In this case, p e g for ZrZn2 indeed 
has a maximum at the Lifshitz-transition point s = pl 
as a function of B. Hereafter we focus on this Lifshitz- 
transition point and perform mean-field calculations by 
using the DOS estimated from fitting to the LDA re- 
sults, 19 which is illustrated in Fig. 3 (b). The fitting 
of the DOS is done by assuming the expansions of the 
DOS given by Eq. (6) and (7). Here we ignore the fine 
structures of the DOS for e > pl corresponding to other 
Lifshitz transitions. Since the qualitative features of the 
MQCP do not depend on the precise location of the para- 
magnetic Fermi level po, we choose po at pl — 2.7[meV] 
to be consistent with the critical value of the magnetiza- 
tion expected from the experiment discussed below. 

Although the MQCP emerges strictly at T = 0, the 
MQCP universality is observable even at finite tem- 
peratures beyond a crossover from the Ising one. The 
crossover is clearly seen in our calculated result for the in- 
dependence of x _1 at a fixed U /W and T. As is described 
above, the TO-dependence of x _1 is in sharp contrast with 
the Ising criticality: Within the lower-field phase, x -1 
shows a convex dependence on <5m and \d\ ~ 1 /dm\ di- 
verges as \m — to^I -1 / 2 (see Fig. 4), whereas x _1 for the 
Ising critical point should have a concave dependence on 
Sm with vanishing \d\ /dm\. Even at finite tempera- 
tures, \d\ _1 / dm\ is largely enhanced away form the crit- 
ical point as is shown in Fig. 4. Here we choose the con- 
trol parameter U corresponding to T c lower than 1.5K. 
As thermal fluctuations scaled by T c become suppressed, 
the enhancement of \d\~ l / dm\ becomes significant for 
lower T c as is shown in Fig. 4 (b). 

In Fig. 4 (c) and (d), the magnetization curve of ZrZii2 
at P=2.1 GPa and T =3.0 K reported in Ref.l (see Fig. 
4 (c)) is compared with the mean-field result calculated 
at T = 3.0 K for the parameter set corresponding to 
T c = 1.5K. Since the first-order transition appears to 
close just near this pressure, we have taken the crossover 
value B* which gives the maximum slope. Because 3.0 K 
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Fig. 4. (a) Inverse susceptibilities x~ 1 as a function of m calcu- 
lated at fixed temperatures T = T c . The thick solid curve repre- 
sents x —1 around the critical point at T = T c = K, namely the 
MQCP. The thin dashed curve represents x -1 at T = T c = 0.15 
K and the thin solid curve represents \~ 1 at T = T c = 1.5 K. (b) 
Absolute values of slope of inverse susceptibility \d\~ 1 /dm\ as a 
function of m. The notations of the curves are the same as (a), (c) 
Measured magnetization curve 2 at P = 2.1 GPa and T = 3.0 K. 
The arrows indicate m* and B* . (d) Logarithmic plots of exper- 
imental and theoretical magnetization curves within the lower- 
field phase; the dashed and solid bold curves represent the ex- 
perimental and theoretical magnetization, respectively. Straight 
solid lines represent 8 = 3/2. The dashed straight lines represent 
the linear dependence; m — m* oc B — B* . 



seems to be still higher than T c in the experimental data, 
we expect that m — m* | oc \B — B* \ when B is too close to 
B* . The experimental data clearly show a crossover from 
\m-m*\ oc \B-B*\ to \m-m*\ oc {B-B*^/ 5 with 6 > 1 
at a very small value of \B - B*\ ~ 0.003 T. The fitting 
of the experimentally observed scaling shows 8 ~ 1.5 for 
\B — B*\ > 0.02 T in remarkable agreement with the 
present prediction in Eq. (9) for MQCP shown as the 
solid curve in Fig. 4 (d). This suggests that the MQCP 
characterized by 6 — 3/2 must be located close to the 
point at (m,B) — (m*,B*), where the fitting suggests 
m* =0.018 jU B and 5*=0.066 T. It would be interesting 
to tune P and B to find the pinpoint of the MQCP. 

An experimental Stoner factor estimated by S = 
x\b=o /2p(Mo) becomes extremely large as S ~ 10 3 , 
which appears to be beyond the applicability of the 
present mean-field treatment. This may originate from 
density fluctuations and may cause a parallel shift of the 
curve as is seen in Fig. 4(d), while the criticality and the 
crossover analyzed in this letter should remain valid. 

The real transition expected all through along the 
QCL instead of the crossover may be detected by a jump 
of the thermopower 20 in this quantum topological meta- 
magnetism, although at low T ^ 0, the QCL becomes a 
sharp crossover. 

We finally discuss the qualitative reliability of the 
present mean-field theory around the MQCP. There are 
two possible limitations of the present theory; an in- 
stability of the homogeneous phase and influences of 
density fluctuations apart from the largely enhanced 



Stoner factor discussed above. Around the MQCP ex- 
pected for ZrZn2, the linear instability of homogeneous 
phases signaled by the negative values of the uniform 
charge compressibility k oc (1 — [/p c ff) _1 , where p e g = 
■J p(^i)p(^i), 13 occurs in a tiny field range given by 
\B-B L \/B L < 5 x 10~ 4 for B < B L and \B-B L \/B L < 
10~ 2 for B > Bl, even if the costs of long-range Coulomb 
energy for phase separations are ignored. Even when 
such a phase separation is expected, the enhancement of 
\dx /dm\ remains observable outside the region of inho- 
mogeneity, whereas the Ising criticality could be masked 
by the phase separation around the critical point located 
at low temperatures. The density fluctuations do not af- 
fect the validity of the mean-field exponents for B > Br, , 
where the expansion Eq. (8) characterized by a ip -theory 
is on the upper critical dimension d c = d + z = 6. On the 
other hand, the upper critical dimension of the expansion 
Eq.(8) for B < Bl may be higher than 6. Therefore, in 
the close proximity of the MQCP, 6 > 3/2 could be ob- 
tained, although it is not seen in the experiment. 

In summary, the metamagnetic transition in three di- 
mensions accompanied by changes in the Fermi-surface 
topology is shown to belong to an unusual universality 
class. The topological and symmetry-breaking characters 
meet at the MQCP, where the inverse susceptibility x" 1 
has a divergent slope as a function of the magnetization 
m, in contrast with the Ising criticality. The experimen- 
tal results for the metamagnetic transition of ZrZn2 are 
consistent with this prediction, implying the relevance of 
the unconventional universality. 
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